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(57) Abstract: A method is disclosed for using a three-dimensional seismic image of a subsurface earth volume to construct a geo- 
logic model specifying the spatially-varying grain size distribution, porosity, and permeability throughout the volume. The method 
applies to earth volumes composed of water-lain clastic sedimentary deposits and involves, in one embodiment, (a) identifying the 
outline forms of geologic bodies in geologic data; (b) using the outline forms of the geologic bodies to determine the spatially-vaiy- 
ing grain size distribution within the bodies, guided by assumptions about the nature and behavior of the paleoflow that deposited 
the bodies; (c) determining rock properties such as, porosity and permeability within the geologic bodies based on grain-size distri- 
bution, mineralogy and burial history information. 
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A METHOD FOR CONSTRUCTING GEOLOGIC MODELS OF 
SUBSURFACE SEDIMENTARY VOLUMES 

CROSS RELATED APPLICATIONS 

[00011 This application claims priority to U.S. Provisional Application No. 
60/609,057, which was filed on September 10, 2004. 

FIELD OF THE INVENTION 

[0002] This invention relates generally to the field of geologic modeling. 
Specifically, the invention is a method for constructing a model of geologic properties 
throughout a subsurfece volume. 

BACKGROUND OF THE INVENTION 

[0003] A geologic model is a computer-based representation of a subsurfece earth 
volume, such as a petroleum reservoir or a depositional basin. In the oil and gas 
industry, geologic models are commonly used in activities such as determining the 
locations of wells, estimating hydrocarbon reserves, or planning reservoir- 
development strategies. Geologic models are commonly inputs to production flow 
simulations that are used to test development and production scenarios in order to 
optimize return on investment A key parameter in flow simulation is the spatial 
distribution of permeability, which together with the properties of the hydrocarbons 
and other fluid found in the subsurfece reservoir, determmes the producibiUty of the 

reservoir. 

[0004] The geologic modeling process can use many different data types, 
including but not limited to rock-property data obtained firom wells and seismic data, 
as well as structural and stratigraphic surfaces that define distinct regions within the 
model. The goal of flie process is to construct a representation of the subsurfece diat 
is realistic, accurate, and consistent with tihe available data types. 
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[0005] Geologic models may take on many different forais. Most commonly, 



of a three-dimensional array of model blocks (cells), to which geologic and/or 
geophysical properties such as lithology, porosity, acoustic impedance, permeability, 

5 or water saturation are assigned (such properties will be referred to collectively herein 
as "rock properties"). The set of cells constitutes the geologic model and represents 
the subsurface earth volume of interest. Dimensions of the cells are commonly chosen 
so that the rock properties within a cell are relatively homogeneous without creating 
an excessive number of cells. The goal of the geologic modeling process is to assign 

10 rock properties to each cell in the geologic model so that the resulting model is an 
accurate representation of the subsurface earth volume of interest. 

[00061 There are two main ways to populate the discretized geologic volume with 
properties-geocellular techniques and object-based modeling. In the geocellular 
approach, geostatistical estimation methods (which may be either deterministic or 

15 probabilistic) are commonly used These methods take into account distance, 
direction, and spatial continuity of the rock property being modeled. Deterministic 
estimation methods typically calculate a minimum-variance estimate of the rock 
property at each block Probabilistic estimation methods develop distributions of the 
rock-property values and produce a suite of geologic models for the rock property 

20 being modeled, witii each model theoretically being equally probable. The spatial 
continuity of a rock property may be captured by a variogram. A variogram is a well- 
known technique for quantifying the variability of a rock property as a function of 
separation distance and direction. Geostatistical models offa: several key advantages 
in that they can utilize a wide range of existing statistical algorithms, readily 

25 accommodate data control points such as wells, and are amenable to rock property 
modeling and optimization using geophysical constraints such as, three-dimensional 
seismic data. U.S. Patent Nos. 5,838,634, 6,381,543 and 6,480,790 cover geocelluar 
modeling methods embodied ia processing flows which include repetitive 
optimization steps to drive the geologic model toward conformance with geologic and 

30 geophysical data, such as well logs, seismic surveys and subsurface fluid production 
and pressure dataT^Mostxomraerciai geoiogrcTn^^ 



descriptive or static geologic models built for petroleum applications are in the form 
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PETREL, GOCAD and STRATAMODEL, contain a ydde spectrum of geostatistical 
tools designed to fulfill the requirements of the reservoir geologist or engineer. 

[0007] The chief drawback to geocellular models is that they generally do not 
closely replicate structures observed in depositional systems such as rivers, deltas and 
deep-water canyons and fens. This consideration is significant in that the iatemal 
stmcture of the depositional system may have a significant effect on reservoir quality 
and continuity. 

[0008] Object-based geologic models treat subsurfece reservoir volumes as 
assemblages of geologic objects such as channels or depositional lobes. U.S. Patent 
No. 6,044,328 discloses one object-based modeling scheme that allows a geologist or 
reservoir engineer to select geologic objects &om an analog library to best match the 
reservoir being modeled. The appropriateness of the analog is judged by the operator 
of the process based on their geologic experience. Most commercial software 
packages, including PETREL, IRAP-RMS and GOCAD implement objects as 
volumetric elements that mimic channels and lobes using sin5)lified elements based 
on user-deformable shapes such as half pipes and ellipses. 



10009] Object-based modeling is most useful where three-dimensional spatial 
information such as that provided by three-dimensional seismic vohimes is lacking or 

I is of low resolution. The simple shapes provided in existing methods can not readily 
capture the complex spatial information seen in modem three-dimensional seismic 
surveys. Furthermore, most of the current techniques emphasize channels and 
channel complexes as the primary reservoir element. While channelized systems are 
significant hydrocarbon reservoir types, the identification or modeling of channels is 

5 incon5)lete without accurately characterizing the potentially porous and permeable 
sands that fill them. The sands themselves are packaged in lobate bodies much like 
sand bars observable in modem rivers or deltas. It is these bodies which form the 
bulk of siUciclastic reservoirs. Geologic models that honor and take advantage of 
these naturally occurring, fimdamental elements should produce more accurate 

} — subsurfece-models:- 
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[0010] Process-based geologic modeling tools attempt to reproduce subsurfece 
stratigraphy by simulating or approximating the physical processes of sediment 
transport and deposition, building sedimentary deposits in chronological order. The 
simulation results can be checked against subsurfece data and the simulation rerun 
5 using new control variables in an iterative process until approximate correspondence 
with subsurfece data is achieved. U.S. Patents Nos. 5,844,799, 6,205,402 and 
6,246,963 describe three such methods which ranploy diffusion or rule-based process 
models to create basin-scale models with limited spatial detail. 

[0011] Process-based models can typically generate realistic-looking simulated 
10 deposits, but they are not commonly used for commercial geologic modeling because 
it is difficult to adjust the model inputs in such a way as to cause the simulated deposit 
to honor subsurfece data constraints. Unlike geocellular or object-based models, 
process-based models typically cannot be efificiently optimized by computet 
algorithms when new static or dynamic data become available. Finally, modem three- 
15 dimensional seismic surveys provide a spatial framework for geologic models which 
is difficult to replicate with process-based models. 

[0012] From the foregoing, it can be seen that there is a need for a method that 
honors the shapes and property distributions of naturally occurring sedimentary 
deposits but can also be easily tied to available seismic and well data. In the above 
20 described method, process-based, object-based and geostatistical approaches may be 
utilized. Preferably, the method may provide an automated optimization process 
capable of being performed by a conq)uter, resulting m a more accurate model of the 
subsurfece earth volume of interest with minimal additional time and effort The 
present invention satisfies tills need. 



wo 2006/031383 



PCT/US200S/029753 



-5- 



SUMMARY OF THE INVENTION 

[00131 A method is disclosed for constructing a geologic model which specifies 
rock properties within a subsur&ce volume based on seismic data: In this 
embodiment, (a) outline forms of at least one geologic body within the subsurface 
5 volume are identified from seismic data; (b) the grain-size distribution at at least one 
point within a geologic body is estimated from the outline form of the geologic body; 
(c) rock properties including porosity and permeability are determined using the 
determined grain size distribution along with mineralogy and burial history 
information. 

10 [0014] A second embodiment of the method for constructing a geologic model 
which specifies rock properties within a subsurface volume based on seismic data is 
also disclosed. In this embodiment, (a) outhne forms of potential geologic bodies are 
extracted from geologic data; (b) the outline forms are tested, selecting only those 
outline forms that likely correspond to sand-prone geologic bodies (c) the size of the 

IS outline forms is adjusted to make them substantially consistent with additional 
constraints on the geologic model such as interval geophysical properties; (d) the 
spatially-varying grain-size distribution is determined within the outline forms; (e) 
rock properties are determined within the geologic bodies using the determined grain- 
size distribution. In addition, the rock properties can be determined using burial 

20 history information along with the grain-size distribution. 

[0015] A third embodiment of the method for constructing a geologic model 
which specifies rock properties within a subsxuface volume based on seismic data is 
also disclosed. The method comprises: (a) identifying the outline forms of at least 
one geologic body from the geologic data; (b) estimating the depositing flow field that 
25 created at least one geologic body from the outline form; (c) estimating the grain size 
distribution within at least part of the geologic body using a relationship between the 
estimated flow field that created the geologic body and the grain-size distribution; (d) 
estimating rock properties of the sedimentary basin within at least part of the of 
sedimentary basin based on grain-size distributioiL 
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BRDEF DESCRIPTION OF THE DRAWINGS 



[0016] Figure 1 illustrates tiiie initial stage of deposition of a sedimentary system 
by a point-source flow, 

5 [0017] Figure 2 illustrates a latra stage of deposition of the same sedimentary 
system shown in Figure 1; 

[0018] Figure 3 illustrates an even later stage of deposition of the same 
sedimentary system shown in Figure 2; 

[0019] Figure 4 is a flow chart of an embodiment of the invention; 

10 [0020] Figure 5 is a flow chart of a second embodiment of the invention; 

[0021] Figure 6 is a three dimensional view of a seed detected seismic cube; 

[0022] Figure 7 illustrates a map view of the three dimensional volume illustrated 
in Figure 6; 

[0023] Figure 8A illustrates a cross sectional and map views of a potential 
IS reservoir interval; 

[0024] Figure 8B illustrates a map view of a potential reservoir interval; 

[0025] Figure 9 illustrates a map view of a single potential reservoir sand body 
which has been populated with estimated grain-size distributions; 

[0026] Figure 10 is a schematic representation of the burial and diagenetic 
20 alteration of a sample block of reservoir sediment with increasing burial time, 
pressure or temperature; 

[0027] Figure 1 1 illustrates a three dimmsional view of a reservoir interval where 
non-reservoir and low quahty reservoir material are added to the geologic model; and 



[0028] 



Figure 12 is a flow chart of a third embodiment of the invention. 
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DETAILED DESCRIPTION 

[0029] In the following detailed description, the invention will be described in 
connection with its preferred embodiment However, to the extent that the following 
description is specific to a particular embodiment or a particular use of the inv^ition, 
5 this is intended to be illustrative only. Accordingly, the invention is not limited to the 
specific embodiment described below, but rather, the invention includes all 
alternatives, modifications, and equivalents falling within the true scope of the 
appended claims. 

[0030] The deposition of clastic sedimentary bodies typically begins with a flow 
10 of sediment-laden water firom a confined channel, such as a river mouth, into an open 
region where the flow decelerates. The point where the sediment-laden flow is able to 
e3q)and and decelerate is typically known as the uilet Initially such flows expand 
fi:eely and deposit sediment as the flow decelerates. Thereafter, as the deposited 
sediment grows in height, it begins to obstruct the flow field Eventually, the deposit 
15 becomes sufificientiy large that the flow is diverted around the deposit. This results in 
a new path and inlet for the flow field to an open region beyond or adjacent to the old 
deposit The deposition process then repeats, and a second body in the system is 
created. In addition, more than one such body may be actively built within the system 
at a time. Overall, the process produces a deposit consisting of stacks of sedimentary 
20 bodies, which is a useful model for the structure of hydrocarbon reservoirs. 

[0031] As will be understood by those skilled in the art, the flow characteristics 
which create such bodies can be analyzed from Navier-Stokes equations and 
conservation principles applied to a sediment laden flow beneath stationary clear 
watCT that can be entrained, such as turbidity currents in the deep ocean. Such 

25 analyses can be extended to flows without vertical entrainment, such as a shallow 
water layer mder air, and these extensions are also within the scope of the present 
invention. The shape and properties of the depositing flow determines the shape and 
internal properties of the deposited geologic body. Furthermore, the flow field is 
related to both the outline form of the deposit and the grain size distribution at all 

30 points within the deposit 
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[0032] Applicants recognized tiiat (a) the outline forms of geologic bodies can be 
extracted from geologic data such as seismic data volumes; (b) the outline forms of 
geologic bodies are related to the depositing flow field that created the bodies; (c) the 
depositing flow field in turn detemaines the spatially-varying grain size distribution 
within the body; and (d) the grain-size distribution either singularly or combined with 
mineralogy and burial history information, determines the rock properties including 
porosity and permeability at all points within the geologic body. The combination of 
these observations leads to a method for constructing a geologic model, specifying 
porosity and pOTneabUily throughout a subsurface volume based on a seismic image 
of the subsur&ce volume. 

[0033] Figures 1, 2, and 3 depict the process of formation of point-source water- 
lain clastic depositional systems. Figure 1 shows the earliest stage of deposition, 
where a sediment-rich flow 10 is emitted from an inlet 12 into an open region 14. As 
die flow expands and decelerates, it loses its capacity to carry sediment, and deposit 
; 16 is formed breath the flow. 

10034] Figure 2 shows a later stage in deposition of the Systran first pictured in 
Figure 1 . In Figure 2, deposit 16 from Figure 1 has become larger and is labeled 16'. 
Deposit 16' is so large that it now diverts the flow 20. Deposit 16' is no longer 
building, but the diverted flow is now creating two new deposits, 22 and 24. 

) [0035] Figure 3 shows a later stage in deposition of the system pictured in Figure 
2. The flow 30, is now confined to chamiels 32 and 34 which transport it past 
completed deposits 16', 22' and 24', and is building deposits 36 and 38. 

[0036] A first embodiment will now be described. With reference to Figure 4, 
this embodiment involves procedures to determine rock properties witiiin a subsur&ce 
5 region using a seismic image of the region. As illustrated in Figure 4, an outline form 
of at least one geologic body within the subsurface region is idaitified (step 401). 
The grain-size distribution in at least one point within the geologic bodies as 
identified in step 401 is determined (step 402). Finally, at least one rock property is 
determined -using-the-grain-sizeLdistributiQn found in Step_4fl2_(stgp. 403), „The rock 
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properties include porosity and permeability. The individual steps will be described 
in greater detail in the following paragraphs. 

[00371 Step 401 is identification of outline forms of geologic bodies and can be 
accomplished by a variety of known methods from geologic data. The geologic data 

5 can include but is not limited to seismic data, outcrop data or studies, core samples, 
reservoir performance data and/or well log data, for example. Seismic data is the 
preferred data type for determining the outline form of a geologic body. The geologic 
data can be obtained using many different prior art acquisition method that would be 
known to persons skilled in the art. This data may be commercially available through 

10 seismic or geophysical companies as well as previously acquired data that may be 
available both in analog or digital form. All such methods of acquiring the geologic 
data are intended to be within the scope of the claims. 

[0038] There are three main methods to obtain an outline form of a geologic body 
from seismic data. The first method is the manual identification of stratigraphic 
15 surfeces in the data which define the boundaries of stratigraphic intervals and form 
the top and bottom boundaries of geologic bodies. Typically, geologists manually 
identify stratigraphic surfaces by recognizing and identifying seismic reflections as 
horizons which correlate to changes in rock properties including stratigraphic 
surfeces. The details of this process are fiuniliar to persons skilled in the art. 

20 [0039] A second method is tiie automated propagation of stratigraphic surfaces, 
which define the boundaries of stratigraphic intervals. This method is typically 
performed with a con^uter extrapolating known or inferred stratigraphic surfeces on 
the seismic data. Examples include seismic coherency or discontinuity programs. 
The details of this process are femiliar to those persons skilled in the art and are 
i^^u^A.>A in nnmrnprMfll annlicfltions such as VOXELGEO. GEOFRAME. 



25 



SEISW0RKS3D, and PETREL. 

[00401 A third method is tiie automatic idaitification of stratigraphically realistic 
geologic bodies as regions having similar seismic attribute values. The preferred 
embodimentfoL steM^l is automatia identification. (^^eolfigljLbodies based pn.tiie 
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methods described in U.S. Patent No. 6,674,689 and co-pending U.S. Patent 
AppUcation Nos. 10/195,582 and 10/654,138. These methods aUow geologic bodies 
to be detected and catalogued at finer resolution than is possible with surfece-based 
interpretation methods. Furthermore, these methods allow researchers to extract 
5 meaningful stratigraphic information &om. volumes over a broad spectrum of spatial 
resolutions using a hierarchical approach that searches these data at multiple scales. 

[0041] U.S. Patent Application No. 10/654, 138 discloses a method for performing 
stratigraphically based seed detection in seismic data that satisfy certain attribute 
criteria that are stratigrsqphically plausible. The criteria honor the layered nature of the 

10 subsurfece so that the resulting seismic objects are stratigraphically reasonable. 
Seismic objects may be extracted that satisfy the iiqmt criteria. Furthermore, this 
method can accurately determine the size and shape of a specific seismic object by 
growing the Object around a chosen seed point in a seismic data volume. The object 
is grown by propagating the object to adjacent traces if the adjacent race conforms to 

IS the selected criteria. 

[0042] U.S. Patent Application No. 10/195,582 discloses a method for performing 
object-based connectivity analysis in a three-dimensional data volume. This method 
provides the ability to predict the connectivity of a seismic object identified by seed 
detection or similar method. In this method, each identified seismic object can be 
20 tracked for changes in its size, spatial position or connections to other seismic objects 
as function of a chosen threshold value or other user contirolled constiaint if two or 
more runs are compared. This allows the operator to pick a tiireshold comiectivity 
value and track a sedimentary body based on that threshold value. 

[00431 U.S. Patent No. 6,674,689 discloses a method for classifying the 
25 morphology of seismic objects. In this metiiod, geometiic statistics are calculated and 
used to classify the morphology of the seismic object Seismically-defined bodies 
that have shapes and internal attributes indicative of good hydrocarbon reservoir 
potential (or other element of interest) can be evaluated and extracted firom the 
background of uninportant seismic objects. 
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[00441 In a preferred embodiment, a sedimentary body is extracted by growing 
the body from a seed point using stratigraphically based seed detection in a seismic 
data volume. The growth of the body satisfies certain attribute criteria designed to 
make the body stratigraphically reasonable. The size and shape of a specific seismic 
5 object can be accurately determined. Next, an object-based connectivity analysis is 
performed on the extracted body. The identified seismic object may be tracked for 
changes in its size, spatial position or connections to other seismic objects as function 
of a chosen threshold value. Finally, llie extracted seismic objects are 
morphologically classified. 

10 [0045] Step 401 can accommodate two-dimensional and three-dimensional data 
volumes of many types, including but not limited to seismic amplitude (near, 6r and 
full stacks), seismic attributes Goop duration, instantaneous amplitude and phase), 
seismic impedance (band-limited and total) as well as derivative products such as 
seismic &cies and discontinuity data. 

15 [0046] Step 402 involves the determmation of the grain-size distribution at any 
point within the geologic bodies. This step can be acconq)lished by a variety of 
techniques. Some of the techniques have additional data requirements, and the 
prefeired embodiment for this step depends on the availability of this auxiliary data as 

■ 

well as on the nature of the geologic bodies. 

20 [0047] One technique is a method for relating the three-dimensional outline of 
geologic bodies to the internal grain-size distribution. This method is described in co- 
pending U.S. Patent Application No. 10/850,018. This is the preferred embodiment if 
fully presCTved (not eroded) bodies can be identified from seismic and should be 
applied if the three-dimensional form of the bodies can be accurately determined. In 

25 this method, the mtemal properties of a water-lain sedimentary body can be 
determined from the three-dimensional shape of the deposit One embodiment 
comprises solving equations rekiting the identified three-dimensional shape of a 
sediment body to the flow field that produced the body and using the flow field to 
solve for tiie sediment body properties in at least one point within the body. 
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[0048] A second technique is a metiiod for relating the two-dimensional plan view 
outline or contour of geologic bodies to flie internal grain-size distribution. This 
method is described in co-pending U.S. Patent Application No. 10/829,600. The 
method should be applied if only the two-dimensional shape of the geologic bodies 

5 can be accurately determined. Using this method, the properties of a water-lain 
sedimentary deposit may be predicted at any location from a contour of constant 
deposit thickness. One embodiment of the method comprises (a) detennining an 
outline of constant deposit thickness in a measured deposit, (b) determining the fluid 
flow properties at the inlet of the measured deposit, and (c) determining a property of 

10 the deposit at any point inside the deposit from modeling the fluid flow 

[0049] A third method is a method for relating the thickness of geologic bodies to 
the internal grain-size distribution. This method is described in co-pending PCT 
Application No. PCTAJS04/05019. If only remnants of bodies are available or if the 
bodies extend beyond the area of the seismic data then this method should be applied 

15 provided a grain-size distribution measurement is available or can be estimated at one 
point within each body. This method predicts the grain size distribution at any 
designated location wifliin a water-lain sedimentary deposit Initially, the vertical 
thickness of the sedimentary deposit at the designated location is determined as well 
as the vertical thickness and grain size distribution at a second location different from 

20 the designated location. Second, a distance parameter corresponding to flie two 
locations is determined. Finally, the distance parameter, the vertical thickness at both 
locations, and the grain size distribution at the second location are used to calculate 
the grain size distribution at the designated location. 

[0050] A fourth technique involves determining the properties of geologic bodies 
25 penetrated by a well and filling a larger geologic region with syntiiesized geologic 
bodies having a similar distribution of propraties to those bodies penetrated at tiie 
well. This method should be applied when the smallest scale geologic bodies are too 
small to be visible in seismic data, provided a core sample or other measurment of 
body thickness and grain size distribution is available. The method for determining 
30 the properties of geologc bodies penetrated by a we ll is d escribed in co-pending PCT 
Application No. PCTAJS04/08238. This method can determine the properties of a 
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water-lain sediment body from a measurement of grain size distribution and deposit 
thickness at one location in the body. In one embodiment, the method comprises (a) 
determining the flow properties at the measured location, (b) extrapolating the flow 
properties back to the inlet through which the depositing flow was emitted, (c) 

5 determining at least one property of the water-lain sediment throughout the sediment 
body by modeling the flow properties, using the extrapolated flow properties at the 
inlet from step (b) as a boundary condition. Once the sedimentary body properties in 
the vicinity of the well are determined, the seismically-defined geologic bodies 
throughout the system can be filled with smaller bodies having an appropriate 

10 statistical distribution of properties as determined from the distribution of properties 
observed at the well. Hiis method is described in U.S. Provisional Patent AppUcation 
No. 60/558,287. 

[0051] A fifUi technique involves simulating the fluid flows which piroduced the 
bodies and adjusting the bovmdary conditions of the flow mitil the simulated deposit 

15 outline most nearly resembles the seismically observed deposit outline. A suitable 
model for this method is disclosed in a co-pending U.S. Provisional Patent 
Application filed concurrently with the present application tifled "a method for 
evaluating sedimentary properties by numerical modeling of sedimentation 
processes." A method for evaluating rock properties by numerical modeling of 

20 sedimentary processes is disclosed in the application. This technique solves fluid 
dynamics equations in a control volume scheme to simulate flow physics and 
deposition. In one embodiment, this numerical modeling technique involves, (a) 
solving a two-dimensional time-dependent mapview system of equations for at least 
flow momentum, flow height, suspended sediment concentration, and entrainment of 

25 overlying water, (b) calculating net sediment deposition at each mapview location 
using the flow properties, and (c) recording the time-variability of the net sediment 
deposition. As computational power increases, these models may have significant 
predictive capability for modeling the rock properties and potential fluid flow 
elements in models of hydrocarbon reservoirs. 

30_ [0052 1 _ The prefraed method for impleme nting fliis fiflh technique is m 
conjunction with hydrodynamics-based gridding (Hydro-Grids) for geologic modeling 
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of the rock properties including grain-size distribution. This method is disclosed in a 
co-pending U.S. Provisional Patent Application No. 60/584,617. In this method, the 
geologic model cell tops correspond to surfaces of constant time as determined during 
the numerical modeling. The rock properties of the basin are then represented in this 
5 hydrodynamics-based gridding structure. 

[0053] Step 403 mvolves the prediction of rode properties within the geologic 
bodies based on grain-size distiibution. This can be acconq)lished by a variety of 
metiiods, including but not limited to the direct calculation of rock properties from 
grain-size distiibution, mineralogy, and burial history. In tiie absence of grain size 

10 data, current geologic modeling methods commonly use cross plots of porosity and 
permeabUily data from core plugs grouped by depositional facies to estimate 
permeability tiiroughout the model. The values and spatial distiibution of porosity are 
usually estimated through geostatistical metiiods linked to depositional fecies analysis 
from two- or three- dimensional seismic data and well control. Persons skilled in the 

15 art of geologic modeling are familiar with these techniques. 

[0054] A preferred method is described in co-pending U.S. Provisional Patent 
Application No. 60/558,284. This metiiod simulates sandstone deposition, 
compaction, and cementation. From the determined grain-size distiibution, 
mineralogy and burial history of the deposited grains, the sedimentation, compaction, 
20 and cementation of tiie graios are simulated. Properties of the sandstone such as 
porosity and penneability may be calculated from the simulation results. The method 
permits multiple mmeralogies to be simulated during the burial history. 

[0055] A second embodiment will now be described. Witii reference to Figure 5, 
tiiis embodiment involves procedures to determine rock properties witiiin a subsurfece 

25 region using a seismic image of tiie region. As illustiated in Figure 5, outline forms 
of geologic bodies within the reservoir interval are identified (step 501). Each outline 
form is evaluated and discarded if it is determined tiiat flie form is unlikely to be a 
sand-prone geologic body (step 502). The outline forms are adjusted in size to make 
them consistent witii additional constiaints on tiie reservoir (step 503). The grain-size 

30 " " distnbufibn arahy point" witB^ tiie geblogfc body fbims ir deteiinihed (step 504). 



wo 2006/031383 



15 



PCT/US2005/029753 



The rock properties are determined using the grain size distribution (step 505). The 
volume within the reservoir interval that is not associated with a geologic body form 
is filled with low permeability material (step 506). 

[00561 Step 501 involves the identification of the outline forms of geologic bodies 
5 from geologic data. The methods described in step 401 of Figure 4 are the same and 
performed in an equivalent mamier in step 501 of Figure 5. A preferred method is to 
automatically extract the outline forms of geologic bodies tcom a three-dimensional 
seismic data volume. 

(00571 Step 502 evaluates each extracted outline to determine whether it is a sand- 
10 . prone geologic body. This is accomplished by seismic attribute methods and through 
geometric analysis of the morphology of the extracted seismic objects. Those 
extracted outlines which are not judged to be sand-prone geologic bodies are 
discarded. 

[00581 Step 503 is to adjust tiie size of the sand-prone geologic bodies to tie 
15 interval geophysical properties detennined from seismic and/or wells. Typically, this 
step involves altering the body outlines by enlarging or shrinking so that the net-to- 
gross of the reservoir interval, as estimated &om seismic or at wells is substantially 
equal to the fraction of vertical distance through the interval contained in body 
outlines. The enlarguig or shrinking process can involve adding or removing a 
20 constant thickness &om the geologic body outline. Alternatively, a constant 
percentage of thickness can be added or removed at each point More advanced 
techniques involve adding or subtracting cells to achieve close matches to sand body 
shapes known from simulations, studies of outcrops, modem depositional 
environments, and high resolution three-dimensional seismic volumes. Persons 
25 skilled in the art will recognize other methods for adjusting body outlines to match 
interval geophysical properties. 

(00591 Step 504 involves the determination of the grain size distribution in at least 
one point inside a geologic body. The methods described in step 402 of Figure 4 are 
the^ama and perfonnfiimjffi jegiiiy^ntimonainstep. 504 o^^^^ 
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[00601 Step 505 involves the prediction of rock properties within the geologic 
body using grain size distribution information from step 504. The preferred 
embodiment is to use grain size distribution and mineralogy along witii burial histoiy 
to simulate sedimentation, compaction, and cementation of the grains. This 
simulation can provide an estimate of porosity, penneability, and other rock 
properties at points where the grain size distribution is known. The methods 
described in step 403 of Figure 4 are the same and performed in an equivalent mamier 
in step 505 of Figure 5. 

[00611 Step 506 is to fill in the space in the model not included in the geologic 
body forms with low peraieabilily material. This filling can be accomplished by 
conventional geostatistical techniques familiar to persons skilled in the art These low 
permeability zones may contain thin sand and silts, but they commonly contain 
significant amounts of detrital clay (very fine-grained particles deposited with the 
sand or silt grains. This clay often fiUs tiie pores between sand grains, reducing 
porosity and permeabUity of the sand. The clay may be concentiated in layers 
("laminated clay"), or fairly evenly distiibuted through the sample ("dispersed clay"'). 
In one embodiment, clustering algoritiuns from Flo-Pac (as described in U.S. 
Provisional Patent Application No. 60/558,284) are used to distiibute clay in our 
simulated rock. The user defines flie amount of clay to be inserted into the rock, size 
and density of tiie clay 'clusters', and the mode of distribution (laminated or 
dispersed). The clay distribution mode, abundance, and cluster density all affect the 
pore structure, and tiierefore ±e flow characteristics, of tiie sediment. Emplacement 
or deposition of clay, growth of cement and physical con5)action aU act to decrease 
porosity and pCTmeability with increasing burial. 

(00621 A hypothetical example will now be discussed. In this example, a seed 
detection of a seismic cube is performed. Figure 6 is a three-dimensional view of a 
seed detected seismic cube 61 showing the individual prospective sand bodies 63 as 
well as the upper 65 and lower 67 bounding horizons of the reservoir interval 69 as 
defined by an interpreter. 
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[0063] Once the bodies are identified through seed detection, they are classified 
based on their morphologies and seismic attributes. Figure 7 illustrates a map view 75 
of the three-dimensional volume illustrated in Figure 6. The sand bodies 63 from 
Figure 6 are classified on the basis of tiieir morphologies and their seismic attributes 
5 into categories that reflect their relative reservoir potential Illustrated in Figure 7 are 
both sand bodies judged likely to be porous and permeable 71 and sand bodies judged 
unlikely to be viable reservoir bodies 73 based on morphology and seismic attributes. 

[0064] In this embodiment, after the morphologies of the bodies are classified, 
object-based connectivity analysis is performed on the bodies in the three-dimraisional 
10 data volume. Figures 8A and 8B respectively illustrate cross-sectional and map views 
of a potential reservoir interval 81 as defined by previous steps of tiiis exan^le. 
Figure 8A is a cross-sectional view dociunenting a potential stacking pattern in a 
hydrocarbon reservoir 83 containing multiple, vertically stacked sand intervals 85 
between the top reservoir horizon 87 and the base reservoir horizon 89. 

15 [0065] Net reservou: to gross interval ratios (net-to-gross) can be defined at each 
column of cells in the section and at each point on the map as shown in Figure 8B. In 
this illustration there are three ranges of net-to-gross illustrated. The ranges are 30 to 
70 perceat net-to-gross 82, five to thirty percent net-to-gross 84, and zero to five 
percent net-to-gross 86. The map gives one possible net-to-gross constraint using 

20 polygons to divide the reservou: interval into domains of relative net-to-gross 
potential. Such a map may be as simple as Figure 8B where broad polygons define 
zones of relative reservoir potential, or it may show much more complex distributions. 
The more complex models of net-to-gross may have spatially-varying patterns derived 
from trace-by-trace interval attribute measurements from seismic and seismically 

25 derived petrophysical volumes. In cases where the assemblage of detected reservoir 
bodies does not total to the appropriate net-to-gross interval value from the map, the 
assemblage may be ac^usted through iterative modification of the areal extent and 
vCTtical tbickness of the individual reservoir bodies. 

[0066] Next, the grain-size distributions and related rock properties of tiie bodies 
30~ 'ie iaedictei "Figwre'9 iUus^ view of a siiigle"p6lefifial'reseivb^^ bo^ 



i-MKjmi:^ COPY 



wo 2006/031383 PCTAJS2005/029753 

-18- 



91 iUustratiiig the flow inlet 93 and the direction of paleoflow 95 that deposited the 
sand body. The sand body has been populated with grain-size distributions based on 
the expected inpxst grain-size and the shape of the body. The grain size distribution 
carried by the flow at the inlet may be estimated from wells or outcrops in the area of 

5 study. The estimated grain-size distributions are generated using one of liie 
algorithms from the suite of shape to grain-size inversion techniques described above 
for steps 402 and 504 of Figures 4 and 5 respectively. In this figure, the grain-size 
contour ranges include coarse 90, medium 92, jBne 94, and very fine 96 sands, but the 
actual output of the shape - grain size methods is complete grain size distributions at 

10 all points in the body. 

[0067] After inversion, the resulting grain size model is taken to the depth, 
pressure and temperature regime ^ropriate for the interval of study. Burial history 
analysis coupled with knowledge of initial grain mineralogies and grain size 
distributions may be used to predict the changes of the sediments &om burial due to 
15 time, pressure and temperature. New simulation techniques such as FloPac (described 
in U.S. Provisional Patent Application No. 60/558,284) permit detailed simulation of 
the sedimentation, compaction, cementation and clay filling, calibrated to thin- 
sections of the sediments. 

[0068] Figure 10 is a schematic representation of the burial and diagenetic 
20 alteration of a sample block 101 of reservoir sedimaits 103 with increasing burial 105 
through time, witii concomitant inweases in pressure and tempaature 107. The initial 
grain pack block 101 is highly porous and permeable witii low density clays 109 and 
no cement The final reservoir rock block 102 has reduced porosity and permeability 
due to mechanical con^action of the sand grains 103 as well as chemical dissolution 
25 and precipitation of pore-filling cements 104 and compacted clays 109. This burial 
history allows detailed estimations of tiie rock properties throughout the history of the 
sedimentary basin. 

[0069] Figure 1 1 gives a three dimension view of tiie last step in tiiis example 
work flow where non-reservoir and low quality reservoir material 111 are added to 
30 ""Be~finji ge61o^c model using caBbra^d to well iand interval: 
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geophysical control. Also shown are the sand bodies 113 as well as the reservoir top 
115 and reservoir base 117. These techniques provide a detailed model of the 
reservoir interval and can be used to estimate both rock properties and reservoir 
production characteristics of the interval In addition, the techniques described above 
5 can be used to predict the properties of the sedimentary basin and not just the 
reservoir interval as described in this example. 

[0070] In addition to the embodiments described above, a third embodiment is 
further described in Figure 12. With reference to Figure 12, this embodiment 
involves procedures to determine rock properties within a subsurfece region using a 

10 seismic image of the region. As illustrated in Figure 12, an outline form of at least 
one geologic body within the subsurface region is identified in step 1201. The 
depositing flow field that created the geologic body is determined from the outline 
form of the body in step 1202. The grain size distribution within the geologic body is 
determined using a relationship between the estimated flow field and the deposited 

15 grain size distribution in step 1203. Then, rock properties wiflrin the geologic body 
are determined based on the estimated grain size distribution in step 1204. 

[0071] Step 1201 is the identification of outline forms of geologic bodies and may 
be accomplished according to a variety of methods, as described in stq) 401. 

[0072] Step 1202 involves determining the flow field that deposited the geologic 
20 body from tiie outline form of the body. In one embodiment, this step is 
accomplished as described in co-pending U.S. Patent Application No. 10/850,018, 
which is hereby incorporated by reference, where equations are given relating the 
three-dimensional shape of a sedimentary body to the flow field which deposited the 
body. In a second embodiment, this step is accomplished in a manner described in 
25 co-pending U.S. Patent Application No. 10/829,600, which is hereby incorporated by 
reference. In that application, the flow properties at an inlet are determined from a 
two-dimensional map view outiine of a geologic body, and tiie flow field that 
deposited the body is determined from the inlet flow properties. In a tiiird 
embodiment, tiiis step is acconq)Ushed by simulating tiie fluid flows which could have 
30 produced tiie body and adjusting flie boundary conditions of tiie flow until tiie 
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simulated deposit outline most nearly resembles the observed outUne from step 1201. 
A suitable method for simulating such flows is described in co-pending International 
Patent Application No. PCTAJS05/21991, which is hereby incorporated by reference. 
In these embodiments, the flow properties specified as part of the flow field typically 
5 include at least the flow velocity, suspended sediment concentration in the flow, and 
grain size distribution of the suspended sediment. 

[0073] Step 1203 involves determining the grain size distribution in the deposit 
ftom the flow properties. This determination is made using empirical relationships 
between the deposited grain size distribution and the flow properties. Appropriate 
10 empirical relationships are described in co-pending International Patent Application 
No. PCT/US05/21991 and U.S. Patent Applications No. 10/850,018 and 10/829,600.), 
which are referenced above. 

[00741 Step 1204 involves the prediction of rock properties wifliin the geologic 
body based on grain size distribution. Methods for making this prediction are 
IS described in step 403. 
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CLAIMS 



We claim: 



5 1. A method for constructing a geologic model of a subsurface volume, 
comprising: 

(a) identifying at least one outline form of at least one geologic body 
within the subsurface volume from geologic data pertaining to a subsurfece volume; 

(b) estimating a grain size distribution witbin at least part of the at least 
10 one geologic body using the at least one outiine form of the at least one geologic 

body, 

(c) estimating rock properties of the subsurfece volume within at least part 
of the sedimentary basin based on the grain size distribution. 

2. The method of claim 1 wherein die geologic data is selected from one of 
15 seismic data, outcrop data, core samples, well logs, and any combination thereof. 

3. The method of claim 1 wherein die at least one outiine form of tiie at least one 
geologic body is identified by a method selected from one of manual identification of 
stratigraphic surfeces, automated propagation of stratigraphic surfaces, automatic 
identification of stratigraphically realistic geologic bodies as regions having similar 

20 seismic attribute values, and any combination tiiereof. 

4. The metiiod of claim 3 wherein tiie at least one outiine form of the at least one 
geologic body is identified by automatic identification of stiatigraphically realistic 
geologic bodies by a particular method selected from one of performing 
stratigraphically based seed detection in seismic data, performing object-based 

25 connectivity analysis in seismic data, classifying morphologies of seismic objects and 
extracting reasonable morphologies, and any combination thereof. 

5. Tte method of claim 1 wherein step (a) fiirther comprises: 
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extracting a seismic body fiom seismic data by growing at least one 
seismic body from a seed point using stratigraphically based seed detection in 
a seismic data set, 

performing an object-based coimectiYity analysis on at least one 
5 extracted body, and 

classifying the at least one extracted body by. morphology and 
extracting the seismic body if the morphology is reasonable. 

6. The method of claim 1 comprising relating the at least one outline form of at 
least one geologic body to the grain size distribution is selected from the group 

10 consisting of relating a three-dimensional outline of at least one geologic body to an 
internal grain size distribution, relating a two-dimensional plan view outline of at least 
one geologic body to the internal grain size distribution, relating thickness of at least 
one geologic body to the internal grain-size distribution, determining properties of a 
water-lain sediment body from a measurement of grain size distribution and deposit 

15 thickness at one location in the at least one geologic body, filling of seismicaUy 
defined geologic bodies with smaller bodies having a distribution of properties 
observed at a well, sunulating fluid flows which produced the at least one geologic 
body while adjusting boundary conditions of the fluid flows untfl a simulated deposit 
oufline resembles the seismically observed deposit outline, and any combination 

20 thereof 

7. The method of claim 1 wherein estimating rock properties of the subsurface 
volume uses at least the grain size distribution in at least one point within the 
subsur&ce volume. 

8. The method of claim 1 wherein step (c) further comprises, 

25 estimating the grain size distribution and mineral conq)osition of grains in 

sandstone, 

simulating sedimentation of grains from the grain size distribution and mineral 
composition of the grains. 
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sunulatmg compaction of the grains, and 

simulating cementation of the grains. 

9. Tte method of claim 1 wherein step (b) fiirflier conq)rises, 

creating upper and lower grid cell boundaries corresponding to depositiorial 
5 time surfaces determined by simulation of sediment deposition. 

creating side grid cell boundaries, and 

modeling the rock properties of the subsurface volume using grid cells 
associated witii the upper grid cell boundaries, lower grid cell boundaries and side 
grid cell boundaries. 

10 10. A method for constructing a geologic model of a subsurface volume, the 
method comprising: 

(a) identifying at least one outline form of a potential geologic body from 
geologic data pertaining to the subsur&ce volume; 

(b) selecting from the identified at least one outline form only the at least. 
1 5 one outline form likely to correspond to sand-prone geologic bodies; 

(c) adjusting the at least one outline form of the potential geologic body to 
make the potential geologic body substantially consistent with additional constraints 
on body properties; 

(d) estimating the grain size distribution within at least a potion of tiie 
20 potential geologic body, guided by the at least one outline form of the potential 

geologic body; 

(e) estimating rock properties within at least a potion of the potential 
geologic body based on the grain size distribution. 
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11. The method of claim 10 further comprising filling at least a portion of a 
volume within a sedimentary deposit that does not correspond to the sand-prone 
geologic hodies with low permeability material. 

12. The method of claim 10 wherein the geologic data is selected from the group 
5 consisting of seismic data, outcrop data, core samples, well logs, and any combination 

thereof 

13. The method of claim 10 wherein the at least one outline form of the potential 
geologic body is identified by a method selected firom the group consisting of manual 
identification of stratigraphic surfaces, automated propagation of stratigraphic 

10 surfeces, automatic identification of stratigraphically reaUstic geologic bodies as 
regions having similar seismic attribute values, and any combination thereof. 

14. The method of claim 10 wherein the at least one outline form of the potential 
geologic body is identified by automatic identification of stratigraphically realistic 
geologic bodies by a particular method selected &om the group consisting of 

15 performing stratigraphically based seed detection in seismic data, performing object- 
based connectivity analysis in seismic data, classifying morphologies of seismic 
objects and extracting reasonable morphologies, and any combination thereof. 

15. The method of claim 10 wherein step (c) is performed by altering the at least 
one outline form so that net-to-gross of at least one reservoir interval, as estimated 

20 from seismic or well data, is substantially equal to the fraction of vertical distance 
ttirough the at least one reservoir interval contained in body outlines. 

16. The method of claun of claim 10 wherein step (d) is selected from the ffovp 
consisting of relating a three-dimensional outline of the at least one geologic body to 
an internal grain size distribution, relating a two-dimensional plan view outline of the 

25 at least one geologic body to the internal grain size distribution, relating thickness of 
the at least one geologic body to the internal grain size distribution, detennining 
properties of a water-lain sediment body from a measurement of the grain size 
distribution and deposit thickness at one location in the sediment body, filling 
seismically defined geologic bodies with smaller bodies having a distribution of 
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properties observed at a well, simulating fluid flows which produced at least one 
geologic body while adjusting boundary conditions of the fluid flows until a simulated 
deposit outline resembles the seismically observed deposit outline, and any 
combioation thereof. 

17. The metbod of claim 10 wherein step (e) further comprises, 

estimating the grain size distribution and mineral conq)osition of grains in 
sandstone, 

simulating sedimentation of grains ftom the grain size distribution and mineral 
composition of the grains, 

simulating compaction of the grains, and 
simulating cementation of the grains. 

18. The method of claim 10 wherein step (d) further comprises, 
creating upper and lower grid cell boundaries corresponding to depositional 

time surfiices determined by simulation of sediment deposition. 

creating side grid cell boundaries; and 

modeling the rock properties of the basin using grid cells associated with the 
upper grid cell boundaries, lower grid cell boundaries and side grid cell boundaries. 

19. A method for constructing a model of a subsurface volume conq)rising: 

(a) identifying an outlme form of at least one geologic body from geologic 
data pertaining to the svibsurface volume; 

Cb) estimating a depositing flow field that created the at least one geologic 
body from the outline form; 

(c) estimating grain size distribution within at least part of the at least one 
geologic body using a relationship between an estimated flow field that created the at 
least one geologic body and die grain size distribution; 
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(d) estimating rock properties of a sedimentary basin within at least part of 
the subsurface volixme based on the grain size distribution. 

20. The method of claim 1 9 fiiriher conqnising using mineralogy and burial 
history information in combination with grain size distribution to estimate the rock 

5 properties. 

21. The method of claim 19 wherein the estimated depositing flow field comprises 
spatial variability of at least one of flow velocity, flow height, suspended sediment 
concentration, suspended sediment grain size distribution and any combination 
thereof. 
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401 



Identify outline forms of geologic bodies 
within a sedimentary basin; 



I 



402 



Fill geologic body forms with spatially-varying 

grain size distributions; 



^403 



Use grain size distributions and burial history 
to calculate porosity & permeability at any 
point within the geologic bodies. 



FIG. 4 



^501 



Identify outline forms of geologic bodies 
within a sedimentary basin; 



I 



/-502 



Analyze outline forms to select those which 
are likely to be sand-prone geologic bodies; 



503 



Adjust geologic body size to fit interval 
geophysical properties; 



I 



504 



Fill geologic body forms with spatially-varying 

grain size distributions; 



I 



505 



Use grain size distributions and burial history 
to calculate porosity & permeability at any 
point within the geologic bodies; 



/I 



506 



Populate the region outside the geologic 
body forms with low permeability material. 



FIG. 5 
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FIG. 7 
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FIG. 8B 
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1201 



Identify outlines of geologic bodies 



1202 



Determine the depositing flow field associated 
with each geologic body from the outline 




1203 



Determine the grain size distribution within each 
geologic body from the depositing flow field 



1204 



Determine rock properties within each geologic 
body from the grain size distribution 



FIG. 12 
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